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lattice calculations where the number of flavors Nf < 4. 
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1. Introduction 

Pseudoscalar flavor-singlet mesons are interesting for a number of reasons, including the rela- 
tion between the v\' and 71 mass difference and the U(l) axial anomaly [|TJ, 0], and the connection 
to the topological susceptibility. Unlike flavor non-singlet states, the propagator of flavor-singlet 
meson contains disconnected correlators, which are particularly challenging to measure precisely 
in lattice simulations. Hence flavor-singlet mesons present a rich proving ground for lattice sim- 
ulations. To this end a number of simulations have been performed on the lattice using quenched 
Wilson [§], N f = 2 Wilson [g |[ |[ |, |], and N f = 2 staggered fermion formulations [g [RJ. 

Simulations utilizing modern dynamical staggered fermion gauge configurations are of partic- 
ular interest. The library of Nf = 2 + 1 staggered fermion gauge configurations from the MILC 
collaboration QI4 [l3| ] contains extremely light up and down quarks (m q /m s < 0.2), simulated at 
fine lattice spacings (a ks 0.1 lfm and a 0.09fm). These configurations, and other ensembles 
currently being generated provide an unprecedented opportunity to measure mixed operators nu- 
merically determine the masses of both the 77 and 17'. Additionally, the sensitivity of disconnected 
correlators to sea quark loops offers a probe of the validity of the fourth-root trick in the staggered 
fermion formulation. 

With this in mind, we make a preliminary presentation of measurements of pseudoscalar 
flavor-singlet propagators on Nf = 2 + 1-flavor improved staggered fermion configurations. As 
the work is in preliminary stages, we concentrate on some theoretical and technical issues in Sec- 
tions || and|3| In Section |] we describe our early results with some discussion of ongoing questions 
and issues, which we hope our continuing work will resolve. 



2. Theoretical Considerations 

The flavor-singlet pseudoscalar meson propagator is expressed as 

G(x',x) = (2>(x')(y 5 ®l)<7i(*0£^.(*)(y 5 <8>l)qj(x)), (2.1) 

where the (75 ® 1) denotes the meson as a pseudoscalar in spinor-space and a singlet in flavor-space. 
One can group the possible contractions into two classes: Nf connected terms with contractions 
connecting fields at x and x', and Nj disconnected terms with the qs contracted with qs at the same 
space time points. So the full propagator is: 

G = N f C-N 2 f D, (2.2) 

where C and D are the connected and disconnected correlators, respectively. The minus sign is due 
to the additional fermion loop in D. The function C is also the connected pion propagator with 
(75 ® 1) operator. In full QCD, both C and G will have leading behavior that decays exponentially: 

C(0~e~ m * f and G(0~e~"V f . ( 2 .3) 

So the ratio of the disconnected to connected contributions to the singlet propagator behaves as 

Njm = N f C(0-G(,) = { j 

V ' N f C(t) N f C{t) 71)1 v ' 
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in full QCD [Q. The correct behavior of the full propagator depends on the sea quarks having the 
same number of flavors and same masses as the valence quarks. If the det 1 / 4 trick introduces some 
unexpected pathology into the action of the sea quarks, we might expect to detect a deviation from 



the form of equation |2.4 For example, in the limit of quenched QCD R should be linear: 

R(t)=A + Bt. (2.5) 
In non-degenerate flavor simulations one should employ obvious generalizations of equations 



2.2| and ^4j. For example, when Nf is split into N q +N S flavors of light and strange quarks we have: 

G — ► (N q C qq +N s C ss ) - (NqD qq + N^D SS + 2N q N s D qs ) 
Np qq +NjD ss + 2N q N s D q 



R 



'qs 



NqC qq -\- N S C SS 



(2.6) 



With staggered fermions, one must rescale the disconnected correlators by 1/4 with respect 
to the connected correlators to take account the two loops of four staggered valence tastes in the 
disconnected diagrams compared to the single valence loop in the connected diagrams^]. In the 
discussion below such rescaling is implicit in D. 

3. Simulation and Measurement 

We have begun measurement of singlet propagators on the following gauge ensembles: 



N f 


P 


L 3 x T 


am 


N ^onfigs 





8.00 


16 3 x 32 


0.02 


104 


2 


7.20 


16 3 x 32 


0.02 


268 


2+1 


6.76 


20 3 x 64 


0.007, 0.05 


422 


2+1 


6.76 


20 3 x 64 


0.01, 0.05 


644 



The 16 3 x 32 configurations are small test lattices generated locally, while the 20 3 x 64 con- 
figurations are part of the library of MILC "coarse lattices" [p^j. 

In principle there are two choices of flavor-singlet pseudoscalar meson operator available. 
These are the (7475 ® 1) and the (75 <g> 1). The former is a three-link operator, with the quark and 
antiquark sources set on opposite corners of the spatial cube, while the latter is a four-link operator, 
with the quark and antiquark situated on opposite corners of the hypercube.. The (7475 <8> 1) has 
a parity partner, namely the scalar (1 <g) 7475), which contributes an oscillating exponential to the 
pseudoscalar propagator. The parity partner of the (75 ® 1), however, has exotic quantum numbers 
J PC = + ~, and contributes nothing to the pseudoscalar propagator. Hence, only the (75 <8> 1) is 
well-suited for looking at quantities such as the ratio in equation |2.4j . We formulate the operators in 
a gauge invariant way, using symmetric covariant shifts to displace the quark and antiquark sources 
to their respective positions on the hypercube. 

We measure the connected diagrams using standard point sources. We measure the the discon- 



nected diagrams with a stochastic volume source method [14]. We define a source field r\ whose 
value at every lattice site is drawn from a gaussian distribution. Then we solve for 



1 



7V src i x,x 4 =ty,y4=t 



t(0, 



(0 



(3.1) 
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where M is the fermionic matrix and A^i is the staggered meson operator that effects the four-link 
shifts and the Kogut-Suskind phasing appropriate to the y r ® 1 meson. We average over a number 
of different source fields v\ , then compute the disconnected correlator: 

D(At) = {0(t)0(t+At)) (3.2) 

In tests on the small (16 3 x 32) lattices, we found that with as few as 40 noise sources our error 
would be dominated by gauge noise even with several hundred configurations, so in subsequent 
runs we used /V src = 40. We also tested Z2 noise and found that it produced noise errors that were 
systematically larger than that obtained with gaussian noise, as shown in Figure [l[ Previously, 
for Wilson fermions Z2 noise was found to be better than Gaussian for the simplest stochastic 
estimators [|T6|]. 




Q-Q Z2 noise 

D— H Gaussian noise 



0.04 0.06 
l/Ns 



Figure 1: Error on the D(t = 1) for Nj ■ = 2 /3 — 7.2 am = 0.20 on 16 3 x 32 lattices for gaussian and Z2 
noise as a function of the inverse number of sources. 



4. Results and Discussion 

Our results to date show clear signals for disconnected and connected correlators. For the 
2+1-flavor ensembles, we formed the ratio: 

AD qq {t) + AD qs {t) + \D ss 

R{t) = 2c M (0+c M (0 • (4 " 1} 

An example of the resulting function is plotted as the dark curve in Figure ||. Results for both 
2 + 1 -flavor ensembles were qualitatively similar. 

It is clear that our correlators as formulated do not form the components of a singlet propagator 
which is positive at all t, as the magnitude of the disconnected correlators exceeds that of the 
connected correlators for t > 5. Correctly normalized correlators showing this behavior might 
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Figure 2: R(t) on j3 = 6.76, am = 0.007,0.05 on 20 3 x 64 lattices for full time series (black), and four 
subsets (colored). Configurations indexed by trajectory number (6 trajectories per configuration). 
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Figure 3: Time series for D qq {t = 10) on /3 = 6.76, am = 0.007,0.05 on 20 3 x 64 lattices. S ame color 
scheme as previous figure. 



illustrate a flaw in the simulated staggered sea. We have suspicions that the normalizations of 
our disconnected operators may instead be incorrect, however. As we compute the disconnected 
and connected correlators by different methods, we must account for different numerical factors in 
each, and reconciling them is non-trivial. Further tests are being made. 

Further uncertainty arises from long autocorrelations in the disconnected correlators. Figure || 
also shows the ratio plotted with four different subsets of the ensemble, each composed of a quarter 
of the time series. The large difference between bins suggests that one needs far more gauge 
configurations to make a statement about the behavior of the D to C ratio. The autocorrelations are 
not visually apparent from inspection of the time series of the disconnected correlators, e.g. Figure 
|3[ but spikes in otherwise small fluctuations seem to be enough to strongly affect D/C. It is also 
possible that with a longer time series R(t) may settle at a smaller asymptote than is apparent now. 
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With these uncertainties, we do not conclude that the pseudoscalar singlet correlators produced 
thus far cast doubt on the staggered formulation, nor are we at this stage able to make statements 
about the lattice masses of the v\ and 77'. It is evident that we will first have to confirm the correct 
normalizations, and then make measurements on extremely long time series ensembles. We are in 
the process of generating 2+ 1-flavor ensembles of ~ 10 4 configurations on the UKQCD's QCDOC 
machine. Additionally, we are investigating different optimization strategies [T5|] such as dilution 
of the stochastic sources which may decrease disconnected correlator noise. 
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